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SUMMARY 


CCHiniN is a rwilw „ w „, in aob.o.tine rot., for U e 
solution Of linear or noa-linoar constrained optimization probers. 
Tho basic optinization al,oritl,e is the Betiod of faasiM.. 
Directions. The user .asc provide a .ain caUin, pro, ram and „ 
estersal routine to e.alaate the oblective and cor.strai._t functions 
ana to provide dradient inforaation. If analytic uradients of the 
objective or constraint functions are not available, tilij 
inforaation is calculated by finite difference. iihiie tea ptojra. 
is intended pri.aril, for efficient solutaon of constrained 
probleas, unconstrained function tininiaation prsnlons ray also no 
solved, and the conjugate direction netted of Fletcher and Eeeves 
is used for this purpose. This manual liscrit.es the use of 
and defines all necessacy paraeeters. Sufficient inforaation is 
provided so that the progras can be used without special knoelody.a 
of optimization techniguas. Sample problems are includes to help 
the user become fasiliar Kith connlB and to make the pro, no 
operational. 
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CONMIS - A FORTRAN PROGRAM FOR 
CONSTRAINED FUNCTION 1 1 N I H I A 1 1 0 N 


SECTION I 
INTF3DOCTION 

.athe.atic.i problem, it is necessary t* detoraiu, 
the aiqiauB oi oaxiaua of a function of several variables, i mte : 
by various linear and nonlinear inequality constraints. it i; 
seldoa possible, in practical applications, to solve these peonies 
directly, an 3 iterative methods are use! to obtain, the tm.ericai 
elution. B» chine-calculation of tais solution, is, of course, 

desirable and the CONJilN proqraa has been .ley eloped to solve a 
variety of ouch problems. 

COMM is a roFTSJll prograu, in subroutine foe*, foe the 
■iai.ization Of a aulti-variaole funztio,. subject to a of 

inequality zonstraints. The general ainiaizution prob.ea is; 
Find values for the set of variables, X(l), to 
tliaiaize OBJ 
Subject to; 

G,J '- LE - !! J=i,:icj» 

*IB(1) .tE.X(I).LE. VUB(I) 1=1, BDV (.lIDn.JT.C 

vhere OBJ is a qeoeral funztion (obloztive funztion) o£ t > 2 
variables, 1(1, , referred to hereafter as lazrsion variables. :; 0 , 
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not be a si.ple analytic function, and . a, ilo any faact 
which can be numerically evaluated. 

G<J) iS the '• 1 “ »* th - inequality constraint. which i , 
also a function of the X (I) . 1C0S is to a na.nar of zoustraiats. 
R(J|. NCON nay be zero. VLDCi uni VU3(I) ^ ^ 

bounds respectively on eariuble X<I). and .ce referred to , s ^ 
constraints. 1«SID2 = 0 indicates that no lower or „ BW: bounU 

Prescribed. I£ «co»-0 and PSIDE-D, tha ob-jecUve function U sail 

to be unconstrained. «» is tha total nu.bar of uociaion 

variables, X(I). 

Co.otr.int t,J) is dBfi „ ed a s actire if CT.«. 8<J , Pr , , 

«d violated, if G (J) . jf. AOS (Cr) , .here constraint thickness. :t. al 
a •■•a specified neqati.e nu.bsr. The nu.ericai sijr ifacanco oi 
cr nay be understood by refer!., to Fi,. 1 , which shows a siuuia 
constr.iat i, . t <0 s.riaM. desiqn sp.co. fonstramt 3 (Jt as 

natheeatically equal to zero alon, a sin, la lino in the iesi-j. 
spaca. Howewer, on . li,ital conpater. the exact „a luo of 3 ,j> - 
can saldoa ba obtained. Therefore, a constraint thickness oi 
2*ABS(CT> is chosen to define the reqion over which r.p, is , SSUI( .j 
to be zero for purposes of opti.i ration, because all s,J, ayst be 
neqatawe, CT is taken as a neqative nuabar for consistency sa that 
•«T GCJ).GI.;t is defined as active ,or violate!) if 

GWI.ST.AflsICT). Khile it eay sees loqical to choose a vet, snail 
»alue (in a.,nitude) for CT (say -1.0E-h), the nature of the 

opti.ization alqoritha used by cosma is such that .ore nuaeticai 

stability can be achieved by ta kiaq ;r=-;. ) or even - a. 2. r is 

used for su.eric.l stability oniy, and when the ooti.izatlon is 
coaplete, one or .ore constraints, ?, (j) , will usually ha .,oai 
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zero, as seen in the examples in GEcrioti VIII. 

It is desirable that 3{J) he normalized so that 
-I.LE.G(J) .LE.1 J= 1 , Nro:« 

In this way the constraint thickness, CT, has tie same numeral. 
sijnificar.ee for ah G (J) . it is not necessary that all 3 (j, u , 
Precisely in this form, and such normalization may not no possible. 
However, it is important that all 3 (J) at least La of the san 
°tder of magnitude. For example, assume that some 3 <J> 
~<ll) **2-X{2) . If X ( 1 ) and X( 2) ace expected to j- of order 1C3 
for the particular problem under consideration, 3 (J, may be scale! 
by dividing by 1C ,'J00 to provide a value for G(j) of oruer one. 

The basic analytic technique used by COKftlR is to minimize obj 
until one or more constraints, 3 (J>, become active. i ha 
minimization process then continues by following the constraint 
boundaries in a direction such that the value of odj continues to 
decrease. Shan a point is reached such that no further decrease i„ 
OBJ can be obtained, the process is terminated. The value of the 
constraint thickness parameter, CX, and the normalization of tne 
constraints, 3(J), have consideraole effect on the numerical 
stability and rate of convergence of the optimization process. 

An example of a constrained nonlinear problem is tl .e 
minimization of the four variable Hoseu-Suzuki function (ref. 1 , ; 

HINIIIIZE OBJ * X ( 1) *♦ 2 - S-X {1> ♦ X(2)»*2 - 5*X(2) + 

2*X{3)**2 - 2 1*X (3) ♦ y (U) **2 ♦ 7*X{4) ♦ s: 

Subject to; 

G{1) = X(1)**2 ♦ X (1| ♦ X (2) *»2 - x ( 2 ) ♦ 

X (3) *»2 ♦ X (3) ♦ X ( 4) --2 - x ( 4) - 


8 . LL. C 


:‘a.J e S 


■ ■■ ' 
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0(2) = X ( 1) * * 2 - X(1) ♦ 2* X ( 2) » * 2 ♦X(3»*-2 

2*X (4) -«2 - X (4) - 10 I 

G(3) = 2* X (1) **2 + 2*X(1) * X(2)*-2 - X (2) + 

X ( 3) +*2 - X (4) - S • , 

• L. _ » ^ 

This: pc»l»l«. h«3 four Secisiot variables , ui t j, lo „ 

constraints, {!.DV = M, NCUN=3). tip loner or upper boui.us VLE(I) ot 
VIID(I) are prescribes so control pari.etar »sip E is spacifiei as 
NSIDE-C to indicate this. It is necessary to provide a sat oi 
initial values for X(i| . and fro. this the constraint optiau. is 
obtained b, C0K„I« and its associated routines. mis probico nrli 
be solve! using COSMIh in SECTION Till. 

The minimization algorithm is based on Zoutenii jk * s method ot 
feasible directions (ref. 2). The alqoritha has bsen modified to 
improve efficiency and numerical stability anti to solve 
optimization problems in which one or more constraints, G(j), are 
initially violated (ref. 1) . While the program is intenaea 
primarily for the efficient solution of constrained unctions, 
unconstrained functions may also be minimized (NCDN=0, and 
NSIDE=0) , and the conjugate direction method of Fletcuar anu ne-vos 
(ref. 4) is used for this purpose. if a function is to 

aaxiaized, this may be achieved by minimizing the negative of the 
function. 

For constrained ainiaization problems, the initial design need 
not be feasible (one or more 3 (J) may be greater than AiiS(CTn , ana 
a feasible solution (if one exists) is obtained vitn a imioal 
increase in the value of the oh-jective function. 

The user Bust supply a main program to call subroutine roKhlb 
along with ar external subroutine to evaluate the objectiva 


function, constraint functions and the analytic gradient or the 
ot jective ar.d currently active or violated constraint 
At any giver. time in the n ininazat ion process, gradient information 
is required only for constraints which are active or violated 
(G(J).GE.CT). Gradients ate calculated by finite difrerenco if 
this information is not directly obtainable, and a subroutine is 
included with COMtilN for this purpose. 

Tha basic program organization is described here, ana 
sufficient information is provided so that the program may oa used 
without special knowledge of optimization techniques. rhe various 
control parameters are described and the required dimensions of ail 
arrays are given so that the user car. limit storage requirements to 
that necessary to solve his particular problems. Sample problems 
are included to aid the user in making the program operational and 
to gain familiarity with its use. 

A summary of all parameters used by conmin and its associated 
routines is given in APPENDIX A for convenient reference. 

APPENDIX B contains a brief diseciption of the subroutines 


associated with CONMIN 
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SECTION ir 

MAKING CUIJMIN OPE9 ATIONAL 

CONMIN utilizes iterative techniques to numerically obtain ti.s 
of • «••«■*! Inaction, so b-) net to m sot u, 

linear and/or non-linear constraints, bocausn of tha natuu oi 
these techniques, the efficiency of tha optimization process (tI ., 
number of times the functions oust ua evaluate!) can oftan i,. 
improvai by the proper choice of control parameters to oca; 
effectively «ith a given problen. For this reason, it is 
particularly desirable that the na. user sol„e several si.ple tao 
to five variable proble.s, experimenting «ith different control 
paraneters. rhe following steps are suggested to help tha user 
•ake CO.ni, operational and to gain the familiarity necessary tor 
its efficient application to a particular problem. 

1. Obtain the deck, including example problems. 

»eai SECTION VIII (EXAnti.ES) of this manual. 

Solve the example problems using COhnlN. and yetiry the 
results by comparison »ith the output listed in this aanual. 

Note that the precise numerical values may differ slightly or. 
different computers. 


2 . 

3. 


**. Read this entire manual carefully. 

5. Devise several too to five variable unconstrained ana 
constrained minimization problems and solve them using zjh-.UK 
If the precise optimum can be determined analytically, compare 
this with the optimums obtained using cosaiN. 

6. Experiment by starting fcon several diffarant initial points 
(different initial x vectors). This is good practice u ail 
optiuizatior. problems, since it improves the chances that the 
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absolute .ini.u. is obtained (instead of a relative * inr.ua, . 
Experiment .Ith various »n»l r tic gradient options 
the same Peebles. .ith aud without calcuiat i,.g 
analytic gradients (see examples 1 and 2 of 5ecti„. , JT;| . 

8. Experiment wih yarious convergence criteria. mm# 
and ITRM. 

•>. Experinent .it,, various constraint thinness parameters, u. 
CT«I», :n and Ctm* to understand the effect of these 
parameters on constrained minimization problems. 

10. Experinent .ith yarious values „ f the push-off factor, t „ Ltk , 
on several constrained sini.ixatiou prohle.s. snail values ot 
THET. .ay be used if constraints. are nearly Unear 

functions of the decision variables, x ( l), and larger yalues 
should be used if one or nor. S ,d, are highly non-linear. 

11 * EXOerl ' ent ^ — lit, options by using „o scaling, automatic 
scaling, or user provided scaling options. 

12. Experinent with yarious conlugate uirectict restart 
pnraeetors, ICNDIE, using oxaepies of unconstrained 
.inlnixation. Pete that if KMI*.,. the steeoest Usctnt 
nethod .ill be used throughout tho opti.izatlon process. 

13. Re-read this entire aanual. 

The default options on the various control para.et.rs hay., 
been chosen as reasonable values for nost optUUation problans. 
The steps listed aboye are i.tenjed to pro.ide the user .ith the 
experience necoss.ry to change those para.eters as regained to 
efficiently solve nee optl.it.tio. problexs of special intenest. 
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SECTION III 
PPOGFUM 3ES ANIZAPI ON 

The overall progra. organization is ohosa aci.aaatic.llv 1: 
Flo. 2. The variables are UitUlUH ir. tna , Jia 
vector x is iteratively changed in subroutine ot„ n , , t< 

associated routines, at.d external routine suit calculates r,„ 
required function values and gradient Information. 

subroutine COM!, is called by the main pcogrw, by too sou 

statement; 

CALL CONSIN (SUB1,0BJ) 

where SOill is the name of the user supplied external subroutine ana 
OBJ is the optimum value of the oblective function upon return from 
cousin. The variables. x,l,. contains in vector x am constraint 
values. S(J>. contained in vector 0 -ill correspond to this option, 
upon return. if adiitional information calculates ir external 
routine SHI is required, the routine should be called agar;, upon 

return from C3NUIH to incurs 1 - 1 .=.*. *.t.- *. 

’ insure that this information correspnd: i3 

the final values of X(I). 

subroutine SUE. is called b, com* and its associate! 

routines by the call statement; 

CALL S0B1 {INFO, OBJ) 

where I.PO is a control parameter defining which information ,ust 
be supplied, and OBJ is the value of the objective function to he 
calculated corresponding to the current decision variables 

contained in I. IBFO win have a value or fro. or.e to tour to 
ideatify shat iaforaation aust be supplied by S031. 

I«ro*1 calculate oblective function value, o 3J . for current 
variables x. 
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I»FO=2 Calculate oblsctiye function valuo. on, ,„ d constraint 
values, G(j), o*1,NCON for current variables, x. 

fa .cu late analytic gradient of objective function 
corresponding to current variables, X. The objective 
function and constraint values already correspond to the 
curront values of X and need not be recalculated. However, 
other intonation obtained in SUB1 when calculating OBJ and 
r ’ (J) not correspond to X and must be calculated again 

here if it is used in gradient, computations. if finite 
difference control nara meter, NFDG, is set to NFPG=1 in the 
main Program this value of INro will never be considered. 

T»F0=4 For current variables, x, determine which constraints are 
active and which are violated (G(J).GE.CT) and how many 
such constraints there are <NAC = Number of active and 
violated constraints). Calculate the analytic gradients of 
the objective function and all active or violated 

constraints. Values of the objective function, OB.T, and 
constraints, G(J), already correspond to the current 
variables, x, and need not be recalculated. \s its the case 
of INFO* 3, all other information used in gradient 

computations must be calculated for the current variables, 

** If finite difference control parameter NPDG, defined in 

the main program, is not zero, this value of info will 
never be considered. 

Note that INFO* J and INF0=4 are considered only if gradient 
information is calculated in the user supplied subroutine, susi. 

With the exception of the external subroutine name, SUB1, the 
objective function value, o«u, and the control parameter, TJPO, all 
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information used in the opti.iaation process is transferred by 
•oanr. of the following labeled oo.mer. blocts. 

COhhoVCfnhl/I MT»T, hBV, IThU. hCOh.hSIDf, ICHDtS. , S c»t,,.D S , 

1 Fi)CH,FDCKm,CT.CTtI»,CTL,CTiaT»,T!IETi.n I 1 I , ! ,A C ,MlC B X1,DI!Ln't', 

? na RF T JK # Ll:<03.],IT-7 ? j,T TER ^ lcsL ^ ) 

C0PM0»/C(n»2/X (vi) ,0F(.vi) (»2) ,isc (ns) ,rc (n») , s (hv.nj) 

c-USM/mm/s (1 » ,r,i (m,.-, 2( „ 2) <C(M ,, #|IS1 

COI^ON/CNMNtt/VLB (N 1) ,VriB(M1) ,SCAL(N5) 

Tho param.»t*rs and arrays used in the minimization process 

1,111 l,0 “ be a * < ‘»* a * f °U««a by the required array di.ensious, hi 
through »o. if a default »alue is listed, this raise uiu be used 
if a zero .alue is transferred fro. the lain program. 
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SECTION IV 

PARAMETERS DEFINED IN .IAIN D R 0 3 n A "I 

mMT Print CO -‘ COJ - eil printing is done or. file S . 

0; Print nothing. 

1; Print initial and final function information. 

2: 1st debut) level. Print all o£ above pin* 

paraiicters. Print function value mi x-voctor at t -,uc, 
iteration. 

3; 2nd. debug level. print all of above plus all constraint 
•alues, nunbers of active or violated constraints, directio: 
vectors, sove para.eters and niscellanoous information, 
constraint paru.eter, DEM, printed unlar this option 

approaches zero as the optUu. oblective is achieve!. 

«; Conplete debug. print all of above plus gradients 0 £ 
objective function, active or violated constraint functions 
and miscellaneous information. 

»DV Humber of decision variables, x<I|, contained in vector X. 

ITH»X Default value a 10. Haxinua number of iterations ii the 
.ini.ization process, if HFDG.Ey.o each iteration requires 
one set of gradient computations (Info = 3 or 4 > ar.U 
approximately three function evaluations (Info = 1 or 2) . if 

NFDG. GT. ° each iteration req lires approximately Him 3 faction 
evaluations (INFO * 1 or 2) . 

" C °* Nu>bf * r of constraint functions, G(j>. neon may be 

Srde constraint parameter. NSIDE*C signifies that tha 
variables X (I) do not have lover or uoper bounds. NSiDa.ll.j 


NSIDE 
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Signifies that all variables X(I| hay,, loner and upper bourne 
defined by VLb<I, and VUB,I, respectively. tf one or „ ru 
variables are not bounded nhile others are. the valuer, of t,... 
loner ar.3 upper bounds on the unbounded variables M4 t be 

taken as very large negative and positive values respectively 
(it*. VLB (I) =-1.02*10, VUD (I) =1. CE* 10) . 

ICKDIB Default value . ®w. Co, jugate direction restart 
para.eter. If the function is currently unconstrained, (aU 
«(J).LT.CT or »CO».»S1DE.O), Fletcaec-3 eeves con -jugate 

direction .ethod will be restarted uith a steepest descent 

direction e.ery ICNDIF iterations. It Iciiai*.! only steepest 
descent will be used. 

"* CW CMtnl »• decision variables win ba 

scaled linearly. 

NSCAL. LT. C ; Scale variables X(I) .by dividing oy SCAL (1) , 
•here vector SCAL is defined by the user. 

NSCAL. EQ. C ; Do not scale the variaales. 

NSCAL. GT. 0 ; Scale the variables every NSCAL iterations, 
variables are noraalized so that scales X (I) = x (I) / jlg « ( i) | . 
Shen using this option, it is desiranle that .scal-iclbis u 

ICNDIB is input as aonzero, and »SCAL=SCV*1 in ICSbit rs input 
as zero. 

RFDG Gradient calculation control parameter. 

srDS.0; All gradient infor.ation is provided by external 

routine SUB1. This infnr»AH«n 

i«is information may do ca iculateu 

analytically, or by finite diffsrei.ee, at the- usefs 

discretion. 

*FDG=1; All gradient information will La calculated t/ tioita 
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difference in COM!*. 3UL1 provides or.lv function velum, 36J 
and G (j) , J= 1, NCON. 

NFDG=2; Gradient of objective function is provided ,,, 
external routine SUM, and dtalients of active ani .routed 
constraints are calculated D, finite difference 1,. con,,;;;. 

This option is desirable if tie qrudiout of tr.o objective 
function is easily obtained in closed for,,, but irudients of 
constraint functions. G(J), are unobtainable. This option ... 

iuprove efficiency if several variables .are Unite: by u„ aL 

or upper bounds. 

FDCB Default value = 0.01. Kot used in BF0 g= 2. Relative cbalje 

i. decision variable x (I) in calculating finite ditterai.ee 

qradieots, For example, FDCH=0.01 corresponds to a finite 

difference step of one percent of the value of the decision 
variable. 

FDCBn Default value . C.C1. sot used if KFdG • 0. niniru. 

obsolete step in finite difference gradient calculations. 
FDCIIK applies to the unsealed variable values. 

Default value = -0.1. Not used if SwQK=«5I D2=0 . 

Constraint thickness paraaeter. If Ci . LS. G ( j) .le.abs(lt) , 
«<J> is defined as active. If G(J) .ST.ABS (CT) , a (J) is said 
to be violated. If G<J).LT.CT, G(J) is not active. :t is 
sequentially reduced in aagnitude during the optimization 
process. If ABS(CT) is very snail. or more constraint, 

■ay be active on one iteration and inactive on the next, nuly 
to becose active again or; a subsequent iteration. This is 
often referred to as 'zigzagging* between constraints, a vid^ 
initial value of the constraint thickness is desirable for 


CT 


Li<jhly nonlinear problems so that when a constraint oecoi (! , 
active it tends to ramain active, thus reuucir. j ta, .iqzaqqi:., 
problem. The default value is usually adequate. 

CTMIN Default value , O.OC4. Not used if NbbN^niDa*: . :U m„u, 
absolute value of -T considered in the optimization process. 
CTH IN may be considered as ‘numerical zero' since it n,ay not 
be meaningful to compare numbers smaller than cihUT. Ih , 
value of CTHIN is chosen to indicate that satisfaction ot a 
constraint within this tolerance is acceptable. The default 
value is usually adequate. 

CTL Default value * -C.01. Not used if :<cON=N3IDt = 0. 

Constraint thickness parameter for linear and oide 

constraints. CTL is smaller in magnitude than Cl because the 
zig-zagging problem is avoided with linear and side 
constraints. The default value is usually adequate. 

CTLRXN Default value = 3. CGI. Mot. used if NCON=NSIDE=0. Minima 
absolute value of CTL considered in the optimization process. 
The default value is usually adequate. 

THETA Default value = I.C. Not used if NC()K = N3IDt.-C. Mean 
value of the push-off factor in the method of feasible 
directions. A larger value of THETA is desirable it tus 
constraints, G (J) , are known to be highly non-linen , ana a 
smaller value may be used if all s (J) are known to be nearly 
linear. The actual value of the push-off factor used in tna 
proqraa is a quadratic function of each G(J) , varying from b.b 
for G (J) =CT to 4.0-THETA for 3 (J) = A63 (CT) . A value of 
THETA* 0. 0 is used in the program for constraints which ar« 
identified by the user to be strictly linear. rasiA is 


call* a 
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PHI 


a •pusl.-otf factor because it pashes the dosi„n ..a, fro, t„ . 
active constrains itto the feasible region. n„. „, £auU 
value is usually adequate. 

Default value - 5.0. hot used if 

Participation coefficient, usel if a desi< ju 13 iafeaaitle (o„ . 
or .ore 8 ( J,. 0 T.»BS(CT„. Phi is 4 .oasure ot to, t „, 

desiq,, win he 'pushed* to.ards the feasible re., lot at,, is, in 
effect, a penalty para.eter. If i,. „ qivoc „ olle , ( __ 

feasible solution cannot be obtained with the default value, 
PHI should be increased, and the proble. run aqiin. j L d 
feasible solution cannot be obtained with PH1=10C, it is 

probable that no feasible solution exists, the default value 
is usually adequate. 

PiCKX 1 »ot used if nsiDOhcosaC. I plus user's best estimate or 
the .ax iau. nu.ber of constraints (including side constraints, 
and VUD (I) ) .hich Bill be active at an, ,i, H „ tilt J 
the airiairation process. MCRX1 . tu.ber of iovs it array A. 
If net, ever exceeds this value, the .ini.ixatiot process 
»lll be teroinated, an error .essaqe .ill bo printed, arc 
control .ill return to tbe aain proqraa. nachxi «iii „, ver 
exceed NDVtl if all constraints 3 (J) and bounds a ,„ 

rn,Z> *" * reasonable value tor u,c;u, , 

the corresponding diaension of array A, is hlb „ , 

.here tbe .ini.u. of .ill only apply for lar ,o probleas a:.d 
i. arbitrary, based on the observation that even for very 
larqe proble.s (o.er a hundred x,i, and several thousand 
it is unco.eon for .any constraints to be ictive at any 
ti.e in the aibi.ir.tie. process (the optisuu solution is 


c_> 
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soldo, -fully const ruined ' for very Ur.,,. „onli„o„r problons, . 
DELFUN Default value . C.C01. Mi.i.u. relative ch„.,e in ls , 
obqective function to indicate convenience. jf 1(; Ir , n 
consecutive iterations, AU5 ( 1 .C-03J (J-1 , /cej , j, , . LI . 6E; FJ s d 
the current desiqr. is feasible (all 1 (J, .IH.iM (cr, , , Ui , 
mini.izatior. process is tar.inated. if tie current desrcn u 
ir. feasible (so.e G (J) . gt. BBS (CT| | , five iterations Jt „ 

required to ter.inate and this situation indicates that a 
feasible desiqr. may not exist. 

D»DPUN Default value = 0.001 ti.es the initial function value. 

Sa.e as DELFOB eicept coaparison is on absolute chanqe in the 

obiective function, BBS (OBJ (J, -OBJ (J- 1, , , instead of relative 
chanqe. 

LI »0&J »ot used if MC3N=KSIDi2=0. Linear obiective function 

identifier. If the objective, 03J, is specifically known to 
be a strictly linear function of the decision variables, X(I), 

set LIN3BJ-1. if OBJ i s a general nonlinear function, set 

LINOBJ=0. 

ItRH Default value = 3. Number of consecutive iterations to 

indicate convergence by relative or absolute changes, D ELF UK 
or DABFUH. 

X(»1) Vector or decision variables, X (I) ,1 =1,«DV. The initial 
X-vector contains the user's best estimate o: the set of 
optimum design variables. 

VLB (81) used only if NSIDB.NE.Q. VLB(I) i 3 the lower ai,owaLie 
value (lower bound) of variable x (1) . if OI .« or «, >re 
variables, X(I), do not have lower bounds, the corresponding 
VLB(I) must be initialized to a very large negative: nu.nl, er 



Page 1 8 


(say - 1 . CF.* 10) . 

VUB(NI) Used only if N5IDE.NE.0. VUB(I) is the maximum allowable 
value (upper bound) of X(I). if Dn( . or iaor , variables. X(I), 
do not have upper bounds, the corresponding vuii(J) swut 
initialized to a very large positive number (say l.dE+10). 

SCAL(N5) Not used if NSCAL=C. Vector of scaling parameters. ; r 
NSCAL.GT.C vector 3 CAL need not be initialize! since SCAL -ill 
be defined in COMBI N and its associated routines. It 
NSCAL.LT.C, vector SCAL is initialized ia the main program, 
and the scaled variables X (I) =x (I) /SCAL (I) . Efficiency of the 
optimization process can sometimes be improved if the 

variables are either normalized or are scaled in such a way 

that the partial derivative of the objective function, dlj, 
with respect to variable X(I) is of the same order of 

magnitude for all x (I). SCAL (I) must bo greater than zero 

because a negative value of SCAL (I) will result in a change of 
sign of Xfl) and possibly yield erroneous optimization 
results. The decision of if, and ho*, the variables should L« 
scaled is highly problem dependent, and soiae experimentation 
is desirable for any given class of problems. 

IS C ( K 8) Hot used if KCOH*0. Lineic constraint identification 
vector. If constraint G (J) is known to be a linear function 
of the decision variables, X (I) , ISC(1) should be initialize! 
to ISC (I) = 1 . if constraint G(J) is nonlinear ISC(I) is 
initialized to ISC(I)=0. Id sntif ication of linear constrain ts 
may improve efficiency of the optimization process and is 
therefore desirable, but is not essential. if ; (J} is u3t 
specifically known to be linear, set ISC(I)=0. 



REPRODUCIBILITY Of THE ORIGINAL PAGE IS POOR. 


13 


SECTION V 

if RAKETL’h S DEFINED IN EXTEE3AL POUTIN’ E JUdl. 

OBJ Value Of objective function f or the current 

variables, X(I) ,1 = 1 , NDV contained ir vector X. Calculate ;, i5 j 
if I NFO= 1 or I NF0=2. 

G(N2) Not used if NCON=NSIDE=Q . Vector containing all 
constraint functions, G (J) , J=1,NC0N for current doci.ro, , 

variables, x. Calculate G(J) , J=1 ,ncon if lH5o=2. 

Analytic qradient of the ob-jective function for the 
current decision variables, X<I). DF(I) contains the partial 
derivative of OBJ with respect to x (I) . Calculate 

DF (I) ,1=1 , NDV if INF0=3 or INF0=4 and if NFDG=0 or MBFG=2. 

Number of active and violate! constraints (G (J) .GE.21) . 
Calculate KAC if INF0=4 and NFDG=d . 

A(N4,N3) Not used if NCON=NSIDF=0. Gradients o* active or violitea 
constraints, for current decision variables, X(l). A(J,lj 
contains the qradient of the Jth active or violated 
constraint, G (J) , with respect to the Ith decision variable, 
X(I) for J* 1 , N AC and 1=1, NDV. Calculate if Info =4 jnu nFjGO. 

IC(N4) Identifies which constraints are active or violatrd. 
IC(J| contains the number of the Jth active or violated 
constraint for J=1,SAC. For example, if c,(10) is the fust 
active or violated constraint (G(J).Lr.cI, j=l,3), set 
IC(1)*1C, Calculate if INF0=4 and NFDG=C. 


If it is convenient to calculate more inf orr,ation tna 


required by the information control paraneto r , Ihfo, thin nay 3 
done. info identifies the minimum amount o£ inforiattj on wh. jh 
necessary at a giver: time in the optimization proems. it i,, 
necessary to determine which bounds (side constraints) V1P(I) , 
VUn(I) are active because this information is determine Ly zo^iK 
Tne required organization of 50 31 is shown m rij. 3 . ., ot 
that if sca!,=G ' «ch , )£ Pi!T . , is lHW11Ciil . 


and can be omitted 
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ITER 


IfCAL (4) 

S(N3) 

61 (M7) 

G2(N2) 
B(M4, N4) 


SECTION VI 

PARAMETERS DEFINED IN CGNMIN AND ASSOCIATED DOUIINLS 

Iteration number. The optimisation process is itarativ 
so that the vector of decision variables at the Et.i 
iteration is defir/ea by X (K> *X (K-1) ♦ALPHA'S (F» , wnoce 1; 
this case K refers to the iteration number ana tu< 
components X(I) are all change:! simultaneously. alpha U 
defined as the move parameter and is printed it the print 

control IPRINT.GE.3. s is the move direction at iteration 
number J. 

Bookkeeping information. NCAL(1) gives the number or 
times external routine SCIB1 was called with INF0=1. 
NCAL (2) gives the number of times INF0=2. :;cal< 3) gives 

the number of times INF0=3 and NCAL<4> gives the number of 
times INFO=4. 

" °’ e direction i„ the NDY-di aansianal optiuUatioc 

S(I) ^ives the rate at which variable x(I) changes with 
respect to ALPHA. 

Not used if »COW=KSIDS=HSCAL=C. Used for temporary- 

storage of constraint values 3 (J) , J=1,ncun and decision 
variables X(I), 1*1 , sdv. 

Not used if NCON*NSIDE*0. Used for temporary storage, or 
constraint values G(J), J*1,MC0H. 

Not used if NCON=NSIDE=C . Used in determining direction 
vector s for constrained minimization problems, array n 
■ay be used for temporary storage i, external routine oJM . 


c (N9) 


Not used 


Used w it .i 


MSI (N6) 


ir ‘ "COH-BStBE-C. 0 S «1 » it „ b x . 

4-.tet.iBiE, direction vector 5 for cor.straUo.i .i.uuatiou 
problems. 0sBd for temporary storage or vector t . 


N'SCAL.Mn.C. routine SOBI. 

Not used if NC3N=USIDE=0. 
deteiminin-j direction vector 
proble®G. Array KS 1 may be 
external routine SUB1. 


Used witn 


array 


B 


m 


S for 


constrained niniaiz ati on 


used for temporary storage i ; . 
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SECTION VII 

PEQ HIRED DIMENSIONS OF A STAYS 


CO M MON/C N h N 2 / X (Nl) ,DF(N1) , i ( N 2 > , ISC(NB) 
COMNON/CNMN3/S (N3) ,G1 (N7) # G2(N2) , C(N9) , 
C0MM0N/CMMN4/VLb(N1) ,VUB(N1) ,SCAL(N5) 


t 1C ( L iH) f A ( N 4 , I, 3 ) 
"1S1 (N6) , h (N4 # N4) 


Dimensions Nl through N7 are minimum dimensions on the 
and are defined by; 

Nl * NDV 

82 * 1 if NCON=NSIDE=0 

* NCUN if WCON.GT.O and N3IDE=0 

* 2*NDV if MCON=0 and NSIDE.GT.Q 

= MCO»+2*NDV if NC3N.GT.C and MSIDE.3S.C 
S3 = NDV+2 

84 * 1 if JJCON® NSID£=0 

NACMX1 if NCOH.GT.O or NSIDE.GT.O 

85 = 1 if NSCAL=0 

= Nl if NSCAL.NE.O 

86 * 1 if NCOS® N3 IDE=0 

= 2*N4 if 5lC08.Gr. 3 or 331 DE. SI . C 

87 « N2 if NSCAL = 0 

= 8AX (N2 # NDV) if NSCAL.NE.O 

88 * 1 if NCON® 0 

* 8C0M if NCO8.GT.0 

89 » 84 if 8FDG.E0.0 

* MAX (N4, NDV) if JIFDG. GT. 0 


a i r a y 


In 


SECTION VIII 
EXAMPLES 

this section several examples are presented, together , lt 

resalts, to provide a setter understanding o: tne pcajI 

organization. In each case the default values are used for co-tco' 

Parameters unless otherwise noted. rho arrav ,H„. 

«ray (lime ns 101.3 iro 

defined in the conon blocks as follows: 

COBBON /CMK2 , X|S 0 ,.«r ( S 0 ».oni 00 ,. ISt( , 009 ,. IC(slwll(S1 5 
™ /csnii/ s l s ,l «s, (Hoot ,12 (hoc) , 33 , ( 1 qj) 

COBHOS /CBBSU/ »LB ( 51 ?) , vuii < 50 ) , SC 4 L ( 50 ) 

The exanples were solved using a CDC 760C co.puter, Ik . 
nuaencal results obtained using other computers „ r Wjir 
slightly froa those obtained here. 


EX4BPLE 1 - cOhSTEaiBED 80SE.-SUSDB! EUdCrlOl.. so 

INFOBHA TION. 


f * r ’ *■ * • “T *"* B r <*% 


consider the .isi.iz.tion problem discussed ir. 3KJI0t I; 
HIBIBI2E OBJ . X(1)., 2 - 5*X ( 1) . M2 ,.. 2 . 5 . 1(2) . 

2*X(3)»«2 - 21»X(3) . X (4) * »2 * 7» x (4) . 50 

Subject to; 

G(1) = X(?)**2 ♦ X( 1 ) ♦ X(21**2 - x (2 ) ♦ 

X(3)**2 ♦ X ( 3) ♦ X (4) *• 2 - X (4 ) - ;« #LE>3 

G(2) = X ( 1) * *2 - X Cl) ♦ 2*X (2) **2 ♦ X(3)~2 ♦ 

2*X(4)»*2 - X { 4) - ic rr . 

• L t • L 


3t 
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G<3) = 2-*(1).. 2 * 2*X(1| * XU)-2 - x (2) ♦ 

*|3)-.2 - *(4, - 5 . ls .j 

this problem is solved baginnin, , U h an initial x-vector of 

* * ( 1 .C, i.c, 1.0, 1.0) 

The opt iaua dasign is known to be 
OBJ = 6 , 0 C C 

and the corresponding X-vector is 

* = (0.C, l.o, 2.0, -I.C) 

print control para.eter of lm«T=2 is „ s0d >a 

gradients are calculated b, finite difference (» ?D 0 = „ . laa 

■axinua nuabar of iterations is taken as itbax=»c to insure not.ai 
taraination. The variables are not scaled, so bsci-s. The 

objective function is nonlinear, so libibj.o. The control 
parameters are defined as: 

IPB1BT*2, BOB-*. IT.dAX.40, BCOB-3, SPUS-I. » ACM ,, 5 1. 

HSIDE=lCNDIR=MSCAL=LISOBJ=ITHH=e. 
FDCB.PDCffN*CT=CTHI3=CT| 1 .CTI.BISi=THeTA = P!I I = D ELFUS.PAtlPulB3. 

BIX constraints are nonlinear so tha linear -on . 

tne linear constraint 

identification vector contains all zeros: 


ISC (J) 


.T= 1 , SCON 


The „i„ prograa and the external routine. EXdPt, ire listed 
i. Pigs. . ,„d 5 respective!,, vith the optimization results it 

Fig. 6. An opti.u. design Of 03J.6.C3J is obtained vith the 
corresponding decision variables; 
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X * (0.0125, G.9°69, 1,9941, -1.C006) 

Bote that the unconstrained minimum of this function „ ay „„ 
found by sotting NCDN =0 in the .min pco,r,.. The undo., strained 

Bininua of OW.-JC.O nay n. found i„ this .ay. and this xs ioft « 

an exercise. 

*n additional problem of interest is to set NCON--2 and, having 
found the optimum subject to these first two constraints only, 
increase SCON to 3 and call C0NdIN again, to obtain th< f if ,al 
optinun design. This is easily done by initially setting NCOH=2 in 
the Min progran, then innediately after returning from cONrtis. aet 
ICOI-3 and call COHRih again. It is not necessary to reinitialize 
the control paraneters. This exercise iomonstra tes the capability 
of COMII to deal *ith initially infeasible dosigns, and suca aa 
option aay ba desirable when minimizing functions for *hich one or 
■ore constraints are difficult or tine-consuming to evaluate, i, 
this ear. the optimization problem may be first solved by ignoring 
constraints ehich are particularly complex. These constraints ate 
then checked to determine if they are violated. if ,. ot , the 
optimization is complete. if one or more such constraints 3 re 
violated, they are added to the sat of constraints, a (J) , =mj 
CORHIR is called again to obtain the final optimum design. Ihis 
approach cannot always be expected to be most africient, but does 
■erit consideration, especially when only moderate constraint 
Violations are expected. 
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RTAUPT.E 2 - CONSTRAINED ROSE!!- SHE [?KI PHNCTION WITH ANALYTIC 

GRADIENTS 

Tb* function .initio i„ FTi-niP 1 is so l,ed hy co , p ,, nn<I 

•11 analytic oradients in closed for.. ,11 control o,ra,eters are 
the MM as before except that no» NFPR*') instead of uppr.,. - he 

uradiont o' the ohiective function vith respect to the decision 
variables is; 


del {OBJ) 



Th© gradients of the constraint functions are; 

( 2*x { 1) ♦ 1 ' 

del(G(1)) « J 2*X(2|-1 

2*X (3) ♦ 1 
V 2*X (4) - 1 


del (C, (2) ) = 



del (G (3) ) * 



The control program is the sa.e as before (Eip. 4 , except no „ 
NEDG=d, and the external routine is naaed EXWP2. Also, the print 
control is now taken as IP*T*T»1 to prowide only su.aary 

inforsation. External routine, EXNP2, is listed in Eig. 7 and the 

optimization results in eia n . . 

m ig. “ t where an optima design of 
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O8J*6.0066 is obtained with 

X = (-0.C072, 1.0022, 2.0053, -0.99C7) 

Tha additional exercises described in example 1 may also De 
solved hare, just as before. 


EXAMPLE 3 - 3-UAR TRUSS. 

As a final example, consider the 3-bar truss shown in Fiq. j, 

Tha structure is subjected to two symmetric, but independent load 
conditions, PI and P2, as shown. The truss is to be designed foc 
siniana weight, sablect to stress limitations only, so that; 

—15 KSI .LE. SIGIJ .LE. 2C ESI 1*1,3 J=1,2 

where SIGIJ is the stress in membar I under load condition j. 

«hil8 this is a very simple structure, it is of particular 
historical significance in the field of automated structure 
design, having been first used by Schnit (ref. 5) to demonstrate 
th** optimally designed structure may not be fully stressed. 

That is, one or more members say not be stressed to their maxima ' 

design stress under any of the applied load conditions. 

The design variables are chosen as the member cross-sect iotal 
areas, A1 and A2, where A3*A1 due to symmetry. Then the objective 
enaction is; 

OBJ * 10*RHO* (2*SQHT (2) *A1*A2) 

where BHO is the Material density (BHO=0. 1) . The stress state is 

defined by; 
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srflll = SIG32 * 20.MSS6T (2 . , 1.» 2 ) / <2. ., 1 . » 2 . s ,„ (2 . , . 4 , . , 

SI021 * SIG22 = 20.-S28T(2. l *»t/(2..»1.»2.S^r ( 2.,*»1.M ) 

SIG31 = HIGH = -20.»»2/(2."«1**2»S2I7(2.) -AI*A1) 

Feaeabering that -15 kst ip cr-ri , _ 

5 KSI * LE * SI * IJ 20 KSI # tiierp ace six 

independent nanlinear constraints <ri =, 

constraints. The compressive stress 

constraint on ae.bet 1 under load condition 1 is ji.on is . 

-SIG11 - 15.0 .LE. o 


or in aoraalized fora; 


SIG11/15.0 -1 .LE. 0 


3iailarly ; 


SIG1 1/20.0 -1 .LE. 0 
•SIG21/15.0 -1 .LE. 0 
SIG21/20.0 -1 .LE. 0 
•SISJ1/15.C -1 .LE. 0 
SIG31/20. V -1 .LE. 0 


Bsc.nsa n.,ati»e .e.ber areas ate not physical!, seaniu,ful, 

Io..r bowds of zero nnst be i.p OS ed on the desi,n sariables. 

«o..»er, notin, that the stress, SIGH, is undefined if »i sqail . 

*«o, loser bounds of 0.001 sill b. prescribed. The upper bound; 

•re taken as 1 . 0 BH 0 to insure that these bounds sill , lste r on 
active. 


The Oblectise function is linear in A1 an„ A2 so the linear 
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oblectiT. function Identifier is token as imonj * i. 

The gradient, o* on., is easily calculated so this .in be done 

analytically, .hlle the qradients of the constraint functions are 

calculated by finite difference. Then nnoo = 2 and the qradient of 
OBJ is defined by; 


del (OBJ) = 


2O.0*S0P'T(2.B)*ph O 


19.C 


The print control .ill be taken as lotTUT „ , the default 

raises are used for all other control para.eters. Then the control 
ear a wet era ere defined as; 

IPBIPT-3, BCOX-6, WSTDB* 1 # WBDG*1, »»CI!X1=51, IT»0B.1*1. 

ITflftXslCSDU^esmtsTTBofsO. 


rDCH=BDCBf!-CT*C7'!IW=rTt=CTLf!TW=TPBTJI»?HI«DBl.POW=0llBFlfX*0. 

*U constrai " ts «e nonlinear so the linear constraint 
identification vector contains all zeros: 

ISC(J) * 0 J * 1,WCO!* 


The lover and upper bounds are defined as 
n,DfT| * 0.9*51 FOB (T) = 1. OB*- 10 

The optiaaa design is known to be 
OB.? * 2.63d 


vhere 
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X * (3.799, C. 408) 

The aain program for this example is listed ir. ?ig. 
the corresponding external routine, EXSPJ, i fl Fi:j . 
optimization results are given in. Fig. 12, where; 

OBJ * 2.639 

and 

x * (C . 789 , C.bC€) 

*ota that only constraint nuabar 2 (the tensile 
constraint in snaber 1 unler load condition 1) is active. 


10 With 
. Th e 


StC £-33 
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^ AT W PROGRAE 

INITI»LTTE DECISION 
VARTABt.ES AND BEPIRE 
CONTROL PARAMETERS 

CALL CONBTN (SgB1,OP,t) 


COWffTN AND 

associated routines 

CALL SUR1 (INFO, OB I) 



SUPPOOTTNE SBB1 (TNFO,OB.T) 

CALCULATE FUNCTION VALUES 
AND GRADIENT TNEOPHATTON 
AS REQUIRED 


FIO. 2 • PROORAfl organisation. 
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SDBROtfTTNT SUB1 (TMPO, OBJ) 

« C 2 !I !10N IPRTNT.nnV, ITBnX,!»CON, ASIDE, TCVDTp wcr»T NPnr 

1 PDCB,PPCHB,CT,CTflTM,CTL r C'T'L*Tv TFF"A p«t «>r»ri n£r™»' NFDG ' 

X'.'.'ll r (ii) , DPtMi) (*2) ,is C (im ,Tc/iai 

«; ^wss&b 



j. 

[cALCOLATB »¥»LTTTC OBA DIPrTop" OBJ I 

j*wn STOBB TB ¥ ECTOR DP. j 

lDP(I) * PBBTTRL DRRT¥TTT¥P OP OBJ 
I jHisra:? TO xm, T « 1,WDV. 


PTO. 1 - O«»0Awt7ATI0B OP EXTBPRAL ROrTTIffP BHB 1 
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c 

c 


c 

c 

c 

c 


c 

c 


c 

10 

c 

20 

c 

c 

c 


EXAMPLES 1 AND 2 OF CONM IN USER * S KANUA 
MINIMIZATION OF CONSTRAINED ROSES-SOZUK 

i C ”2” / C N ” N V 1 p aiMT * N EV , IT MA X , N CO N , N 3 
f DCH, FDC HM, CT ,CTMIN,CTL, CTLMI N , THET A P 
2 DA3FUN,LINODJ,ITBM, ITER , NCAL f 
COMMON /CNMN 2/ X ( 50) , DP ( SO) , 3 { 1 1 00 ) , I SC 
COMMON /CN-HN4/ VLB (5C ) , VU3 (50 ) SCALfDCl 
DEFINE NAMES OF EXTERNAL ROUTINE^ SU61 
EXTaPNAL HXMP1,2XMP2 
INITIALIZE CONTROL PARAMETERS. 


I.. 

I F tj 3 
IDE, I 
HI ,'A A 


cr: > 

C’i DI 
C, N/i 


k , NSC A 

cexi,d 


^ f )W" D , J , 

ELPJI. , 


( 1 0 1 j ) , I c ( j 1 ) , A ( 5 1 , 6 2 ) 


FOR EXAMPLE 2, 

IPRINT= 1 

IPSINT=2 

NDV=4 

IT«AX=40 

HCON=3 

F33 EXAMiLE 2, 

HFDG=0 

N FDG= 1 

NACHX1=51 

NSIDS=0 

IC»DIR=0 

NSCAL=0 

LINOB J=3 

ITRH=0 

FDCH=0. 

FDCHM=0. 

CT=0. 

CTIJIN=0 . 


-hangs the following 


STATEMl 


CHANGE THE FOLLOWING STATEMENT 


TO, 


CTL=0. 
CTLMIN=0. 
THBTA=C . 
PHI=0. 


DELFUN=0. 


DABFUN=0. 

INITIALIZE CONSTRAINT IDENTIFICATION VECTOR 
DO 10 J=1,NC0N 
ISC (I) =0 

INITIALIZE X-VECT9R. 

DO 20 1=1, ND? 

X CH *1 - 


SOLVE OPTIHIZATIOH. 

FOR EXAHPLE 2, CHANGE THE 
CALL CONMIB {EXMP2,OB.J) 
CALL CONM IN (EXRP 1 ,OBJ> 
STOP 
END 


following 


STATEMENT 


ISC. 




FIG. 4 - MAIN PROGRAM FOR EXAMPLES 1 AND 2. 


n n n o 
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S UBROUTI Nt EXMP1 (INFO, OBJ) 

COMMON /CNHN2/ X (50) # D?(53) ,3(113^1 T . . , 

SOUTINE TO PROVIDE FUNCTION AND CONST! AIN"’’ n 

0f c o«s t «»:» E D tt '“ 

p5ScTIO!iVAUJE?' ,, ' ,!NU3Eii ' :! ‘' ,,BUAI " 

* x * a » *2 .-x ,3, .. 2 -^ . (3 , * 

IF (INFO. Ey. 1) RETURN 
CONSTRAINT VALUES. 

§ *1* J u :: 2:5 « 1 »:*»?»** fr* « 2 » *«««■* 2 -* « 3» .x,*, .. 

END 


2-21. «X (3) ♦ 


PIG. 5 


EXTERNAL ROUTINE, EXHP1, FOR EXAMPLE 1. 




=V 


C 0 N M I N » 

FORTRAN PROGRAM FOR * 

CONSTRAINED FUMCTlON MINIMIZATION # 

* n *SA/AMES RESEARCH CENTER. MOFFETT FIELD. CALIF. * 


CONSTRAINED FUNCTION MINIMIZATION 
CONTROL PARAMETERS 


IPRINT NOV 
* 4 


NACMX1 

SI 


LIN08J 

0 


FDCH 

1.00000E-02 

CT 

-l.OOOOOE-OI 

THETA 

1*00000E*0Q 


IT MAX NCON 

♦0 3 

ITRM 

3 

FDCHM 

1.00000E-02 

CTMIN 

A.OOCOCc-OJ 

PHI 

2.SOOOOE*9I 


INITIAL FUNCTION INFORMATION 
OBJ ■ 3.100000E«-01 


X-VECTOR 
1 • 0000E*00 


NSIOE 

0 


ICNDIR NSCAL 


NFOG 

1 


CTL 

-1.00000E-02 

OELFUN 

i.OOOOOE-04 


ctlmin 

1 • OOoOJt-03 

DAB FUN 
3.10000E-02 


l.OOOOE.OO l.0000E*00 1.0000£*00 


CONSTRAINT values 

•4*0000£*00 *6«0000E*00 -1.0000t*00 


THERE ARE 0 ACTIVE CONSTRAINTS AND 0 violaTEO CONSTRAINTS 


6 - OPT If* T7 ATICN" GUILTS FOP FXA'JPLF 1 
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ITER » 1 

X-VECTOR 
1 • 0436E* 00 


OBJ * 2.84843E.01 

1*0436E*00 1«2479£*00 


ITER * 2 

X-VECTOR 

-6.5496E-01 


OHJ s 1.22043E*01 
1 • 0325E.00 2.3S72t*00 


ITER « 3 

X-VECTOR 

2.2440E-01 


0«J * 0.J7b29£*OO 

V.R26BE-01 2«034SE«00 


ITER * 4 

X-VECTOR 

-3.4392E-01 


OBJ • 6.942O3£«00 

1.Q943E+00 2.1498E*00 


ITER • 5 

X-VECTOR 

-6.7S66E-02 


OBJ » to # 3?7U8E.G0 

1*O136E*0O 2.0734£>00 


ITER « 6 

X-VECTOR 

-9.4S81E-02 


°*J • 6.17226E*00 

*# 9 247E-0l 2.0400000 


ITER * 7 OBJ * 6.07060E*00 

X-VECTOR 

7.4640E-02 9.8928 001 1,94780 00 


FT(J. 


8.6847001 


I.3804E-01 


-3.1841t.01 


-8.0388E-01 


6.1323E-01 


9.6346E-0I 


1 .0b62E*00 


6 - OP7IHI7.XTICX RESULTS 


FOR EX**PLE 1 - CC9T. 


ITER ■ 


8 


OBJ * 


6.0218*E*00 


X-VECTOR 

-1.7653E-92 I • 0038E*00 2.0139E*00 -9.f523E-01 


ITER * 9 

X-VECTOR 

1.2500E-02 


OBJ a 6« 01829E*U0 
9.9690E-01 1»9943E*00 -1.0006E*00 


ITER a 10 


ITER « 11 


OBJ • 6.01829E«00 NO CHANGE IN OBJ 

OBJ * 6.01B29E+U0 NO CHARGE In OBJ 


Fl<i. 6 - OPTIRT7 A7T0S RESULTS FOP FXA1PL? 1 - CCH 


final OPTIMIZATION INFORMATION 
OBJ a 6.0I8287E*00 
X-VFCTOR 

3.2*001-0* 9.969OE-0! K99.3t.00 -1.0006E.00 

CONSTRAINT VALUES 

-1.7146E-02 • l . 044SE ♦ 00 -2.8422E-14 

ACTIVE CONSTRAINT NUMBERS 

NO. GREATER THAN 3 INDICATES SIDE CONSTRAINT 


termination criterion 

ABSlOBJii>-O0J<I-l>> LESS Than OABFUN FOR 3 


NUMBER OF ITERATIONS a n 
OBJECTIVE function was evaluated 
constraint functions were evaluated 


68 TIMES 
64 TIMES 


FTf;. 6 - OPTIMIZATION RSS'JLI S FOP FX AMPLE 1 - 


P VT** 4 1 


I TERATIONS 


CC NT. 


Ptpj t*2 


c 

c 

c 

c 


c 


10 

c 


c 


20 


30 


SUBROUTINE EXKP2 (INFO, OBJ) 

C J“ KQN /CNMN1 / IPRINT, NDV. IT *AX -j^r k i-' - • - - r - . . 

I PDCa.FQCHH.rT.rTMT^ rZ, ^',35CAL, 


1 pnr*. c^-,« ” , “ ' io.jl, if;;, ji 

1 r DLa, Y o .. .1^ T , cTP. I N -CTL . Z1 L Ml K : 

2 DABFfTN. T. .'NAPT— TO, 


MX 1 , D 


» <-* - <- 


) , I (SI J , ,i ( , 1, 


2 DABFUN,LiNOB.T f ITi<f1 / ITER*NCAL(4) # 
cof.noN /cn a t; 2 / X(50) ,df( 53) ,3(iiooi i i'/i 

ROUTINE TO PROVIDE FUNCTION AND CONSTRAIN" v - - ' • - 

“F cossrMiKSD «03ES?sari”‘ hA^'" '' 

EXAMPLE 2 OF CONMIN USER»S MANUAL “ 

IF (INF0.GT.2) GO TO 10 
FUNCTION VALUE. 

i??urL { - ) r 2 ' 5 -‘ X(1 > * X(2) **^-5.*X{2)+2.*;< O) 2-21 *>■(!,♦ 

1X(4) **2*7.*X(4) *50. 1 * K ' j} * 

IF (TNF0.2Q.1) PETUPN 
CONSTRAINT VALUES. 

r i ?i ~y / 1 ! I* " + X ^ +X - * 2_X * 2) +x * J) * ( 3 ) + \ ( 4 ) ^ - w _ 

C/3 \-l Jt, ‘it'll V 2 ** X (2) ** 2 * X < 3 > **2*2. “X (4) «*2-V(4)-j- 

CONTINUE 

GHADIENT OF OBJECTIVE FUNCTION 
DF(1) *1.*x (1) -5. 

DF(2)=2.*X{2) -5. 

DF(3) ='!.*X (3) -21. 

I>F (4) = 2. * X (4) *7. 

IF (INTO. E^j. 3) RETURN 

GRADIENTS OF ACTIVE AND VIOLATED CONSTRAINTS. 

IF (G(I).LE.CT) GO TO 20 
N AC= 1 


IC(1) = 1 

A{1,1) *2.*X(1) M. 

A(1 ,2) =2.*X{2) -1. 

A (1, 3) *2.*X (3) +1. 

A (1,4) - 2,*X(4) -1. 

IF (G (2) . IT. CT) GO TO 3C 
NAC=N AC+ 1 

IF (KAC. EQ.NACMX1) RETURN 
I C (MAC) -2 

A ( S AC, 1 ) = 2.*X(1) -1. 

A (U AC, 2) =4.*X (2) 

A (NAC, 3) =2. *X (3) 

A (NAC, 4) =4. * X (4) -1. 

IF (G ( 3) . LT. CT) RETURN 
NAC=MAC+ 1 

IF (NAC.E^.NAC^XI) RETURN 
IC (NAC) = 3 

A(NAC, 1) =4. *X { 1) *2. 

A (NAC ,2) =2.*X (2) -1. 

A (NAC, 3) ~2.*X(3) 

A (NAC , 4> »-1. 

RETURN 

END 


FIG. 7 - EXTERNAL ROUTINE, EX15P2, FoF EXAM PL.’ 


> 
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C 0 N M I N 

fortkan program for 
CONSTRAINED FUNCTION MINIMIZATION 
NASA/AMES RESEARCH CENTtH, MOFFETT FIELD* CALIF. 


INITIAL FUNCTION INFORMATION 
OBJ ■ 3.1 Q000OE*0 1 

X-VECTOR 

I • 0000E*00 1 .OOOOE+OO 1.0000t+00 

CONSTRAINT VALUES 

-+.OO0OE«0O “6.0000F.*Q0 -I.0000t*00 

THERE ARE 0 ACTIVE CONSTRAINTS ANO 


l.OOOOE+OO 


0 VIOLATED CONSTRAINTS 


rTG. 8 - OpTI^T/^TTC'I TSS'JL'I S FOP FXA1PLE 2. 


final OPTIMIZATION information 
OBJ * 6 , 006646E+00 


je 44 


X-VECTOR 

-7.17J0E-03 1.002?k>00 2. OOBJt^OO -9.9073E-01 

CONSTRAINT VALUES 

-S.0054E-03 -1.00B9E+00 -8.S265E-14 

ACTIVE CONSTRAINT NUMBERS 

NO. GREATER THAN 3 INDICATES SIDE CONSTRAINT 


termination CRITERION 

abs (obj ( i ) -08 j ( i- i j ) less Than oabfun for 3 


NUMBER OF ITERATIONS * 13 

OBJECTIVE FUNCTION NAS EVALUATED 
CONSTRAINT FUNCTIONS WERE EVALUATED 
GRADIENT OF OBJECTIVE WAS CALCULATED 
GRADIENTS OF CONSTRAINTS WERE CALCULATED 


34 TIMES 
34 TIMES 
Yd TIMES 
12 TIMES 


I-. 


I TERATIUNS 


* “OPTIMIZATION FESIILT5 FOR fiXA-PlF 2 - CCNT. 
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C 

C 


c 

c 


c 

10 

c 


20 

c 

30 

c 


SXAMPLE 3 OF CONMIN USER'S MANUAL. 

MINIMIZATION OF WEIGHT OF THE 3-3 AR tr|i« 

1 C 2Sc3%n cTli IP ««r.NDV,ITaAX # NCO^NSlt;;,ICN 

o A?ni ,#FDl Lilf1 ' C - T ' C TMlN,CTL,C7lMIN, THETA, PHI NAC 
2 DABFUN # LlNOBJ # lrRM # ITER,NCAL ( 4 ) # ' 

ossss: sss 

INITIALIZE CONTROL PARAMETERS. 

IPR INT= 3 


013, NSCAL, KH. 
n'ao.ix i , del;- j:, , 

retail ,A (5 1,?;^) 


ND V=2 
NCON=6 


NSIDE= 1 

S FDG=2 

NACMX1=S1 

LINO 3 J= 1 

ITMAX=0 

ICNDIR=0 

NSCAL=0 

IT8B=0 

FDCH=C . 

FDCHM=0. 

CT=0 . 


CTB IN=C . 

CTL=Q. 

CTLMIN=0. 

THETA=0. 

PHI=0» 


DELFUN=Q. 

DABPUN=0. 

0S I T5 A iif%coif STBaZNr ID “ TIrIC *"*™ ViCTDK, 
ISC(I)=0 

E „ L3WEF AND UPPE8 BOUNDING VECTORS. 

VLB (I) = C0.'<)01 
VUB (I) =1. CE*10 
INITIALIZE X-VECTOR. 

DO 30 1=1, NDV 

XII) =1. 


SOLVE OPTIMIZATION. 
CALL COKNIN (EXNP3,OBJ) 
STOP 
END 


ISC. 


FIG. 1C - MAIN PROGRAM FOR EXAMPLE 3 


a n n 
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C 


c 


10 

c 


SUBROUTINE EXMP3 (INFO, OBJ) 

COMflON /CNMN2/ X(50) , DF (53) , G M 1 03) t^,i 
RUUTIN’E to provide function and gradient 
3-JAR TRUSS. A 

EXAMPLE 3 01 CONMIN USER'S MANUAL 
RHO=0. 1 


c r : ) , r c < 5 1 > ,a ( 5 1 „ S 2 i 

INFORMATION FOR 


10 


A 1=X (1) 

A2=X (2) 

IF (INFO. GT. 2) GO TO 
OBJECTIVE FUNCTION. 

OBJ= 1C . * RHO* (2. *S 23T (2. ) * A 1 + A2) 
IF (INFO. EQ . 1 ) RETURN 
CONSTRAINT VALUES. 

DENOM = 2. *A1*A2+SQRT (2.) *A1*A1 

= (SQFT(2.)*A1 + A2) /DENOM 

SIG21=20. *SQRT (2. ) *A1/DEN0M 
SIG31=-2C.*A2/DLNOfi 
0(1) =-SIG 11/15.-1. 

G (2) = SIG11/20.-1, 

G (3) =-SIG 2 1/1 5. - 1 . 

G (4)=SIG21/2C.-1. 

G (5) —SIG31/20.-1. 

G (6)=SIG31/2C.-1. 

RETURN 

CONTINUE 


GRADIENT OF OBJECTIVE FUNCTION. 
DF(1)=20.*SQRT(2.) *RHO 
DF(2)=10.*RHO 
return 

END 


FIG. 11 - EXTEPNAL ROUTINE, EXMP3, 


FOR EXAM RLE 3. 
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C 0 N M I N « 

FORTRAN PROGRAM FOR « 

CONSTRAINED function minimization * 

NASA/AMES RESEARCH CENTER. MOFFETT FIELD, CALIF. * 


CONSTRAINED FUNCTION MINIMIZATION 
CONTROL PARAMETERS 


IPRINT 

3 

NOV 

2 

itmax ncon 
10 6 

NSIDE icnoir 
1 3 

NSCAL NFDG 

0 2 

NACMXl 

26 

LINOBJ 

1 

ITRM 

3 



FDCH 

1.00000E-02 

FDCHM 

i-OOOOOE-02 



CT 

•I • OOOOOE-O 1 

CTMIN 

*.00000003 

CTL 

-I.OOOOOE-02 

CTLMIN 

1.00000E-0J 

THETA 

1 • OOOOOE.OO 

PHI 

2.50000E*01 

delfun 

1.00000E-04 

D ABF UN 
3.82843003 

LONER BOUNDS ON 
1 • 00O0E-03 I 

DECISION VARIABLES 
•OOPOE-03 



UPPER 0OUNOS ON DECISION VARIABLES 
1 .0OOOE+ 1 0 1 .00P0E+10 


INITIAL FUNCTION INFORMATION 
OBJ ■ 3.828427E*00 

X-VECTOR 

l.OOOOE+OO 1.0000E+00 
CONSTRAINT VALUES 

-I.9428E*00 •2.9289E-01 -i.5523E*00 

THERE ARE 0 ACTIVE CONSTRAINTS AND 


-5.8579E-01 -6.O948E-01 - 

0 VIOLATFD CONSTRAINTS 


FI'T. 12 - OPTIMIZATIC?! RESULTS POP EXAMPLE 3. 


.2929000 
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begin iteration number 1 

S-VECTOR 

-2oB2 64E>00 "1 • OOOOE+OO 

SLOPE - -9.OOOOOE.OC ALPHA . 1.16J8SE-01 

OBU * 2.798966E^00 

X-VECTOR 

6,76*7E-01 B.B562E-01 


BEGIN ITERATION NUMBER 2 
CONSTRAINT VALUES 

-?.33XU.0(I -1.6S82E.03 -l.6912t.00 -4.81S8E-01 -3.6011E-0; 

P 9!67I1E-01 CTORS * THETA<I, » J*1*NAC 


ACTIVE CONSTRAINT NUMBERS 

NO. GREATER THAN t INDICATES SIOE CONSTRAINT 


CONSTRAINT PARAMETER* BETA . 8.16624E-04 

s-vector 

2.2667E-0? -9.039BE-02 

SLOPE « -2.62B7IE-02 ALPHA « 6.44333E*0O 

OBJ * 2,629589E*00 

X-VECTOR 

®*2252E-01 3,.0315E-01 


-1.4799E»Q0 


begin iteration NUMBER 3 
CONSTRAINT VALUES 

-2.3633E.00 T.6916E-03 -2.06S6E.00 -2.of.79E-0l -T.2229E-01 

p y s ";"" FACTORS, TMETA(I), i.i.nac 
I • 1554E»00 


-1.2U83E*00 


FIG 


12 - 


OPTIMISATION 


SFS'JtTS POP EYAMPLP 1 - 


CCK’T. 


Y'f tf ttlfiMi l : 


Paqf? SO 


I 

-? 


NO, 


OF ACT I VF CONSTRAINTS a 1 OF WHICH 
0 ACTIVE OR VIOLATED CONSTRAINTS ARE SIDE 


0 ARE VIOLATED 
CONSTRAINTS 


ACTIVE CONSTRAINT NUMBERS 

NO, GREATER THAN 6 INDICATES SIDE CONSTRAINT 


2 


CONSTRAINT PARAMETER# BETA * 2.S3279E-05 

s-vfctor 

-1.0595E-02 2.7380E-02 

SLOPE * -2.58609E-03 ALPHA * 1.05049E-02 

OBJ « 2,629562E*O0 

x-vector 

*«2241E-01 3.0344E-01 


BEGIN ITERATION NUMBER 4 
CONSTRAINT VALUES 

•?,3433E#0O 7.4795E-03 -2.0654E«00 -2.0098E-01 

PUSH-OFF FACTORS. THETA (I). I«I,NAC 
1 j4851E#00 


•7.2205E-0I 


NO. OF ACTIVE CONSTRAINTS « 1 0 F WHICH 0 ARE Vim ATEr. 

0 ACTIVE OR VIOLATED CONSTRAINTS ARE SIDE CONSTRAINTS 


ACTIVE CONSTRAINT NUMBERS 

NO. GREATER THAN b INDICATES SIDE CONSTRAINT 

CONSTRAINT PARAMETER# BETA * 1.63864E-0S 

S-VECTOR 

-9.0708L-03 2.3720E-02 

SLOPE ■ -I.9360UE-03 ALPHA « 4.2548&E*00 

OBJ a 2.621324E*00 

X-VECTOR 

7.8381E-01 4.0436E— 01 


FIG. 12 - OPTIMIZATION RESULTS FOR EXAMPLE 3 - CONT. 


•1 .2085E+00 



BEGIN ITERATION NUMBER 5 
CONSTRAINT VALUES 

-2.3423C00 6.7277E-03 -1.9835E*00 -2.6236E-01 -6.4122E-01 

PUSH-OFF FACTORS* THETA <I>. 1*1, NAC 
1.4320E*00 

NO. OF ACTIVE CONSTRAINTS * 1 OF MHICH 0 ARE VIOLATED 

0 ACTIVE OR VIOLATED CONSTRAINTS ARE SIDE CONSTRAINTS 


ACTIVE CONSTRAINT NUMBERS 

NO. GREATER THAN 6 INDICATES SIDE CONSTRAINT 
CONSTRAINT PARAMETER* BETA * 6.46576E-11 


OBJ * 2.62) 324E*00 NO CHANGE ON OBJ 


BEGIN ITERATION NUMBER 6 
CONSTRAINT VALUES 

-2.3423E*00 6.7277E-03 -1.9833E«00 -2.6236E-01 -6.41 2.2E-01 

PUSH-OFF FACTORS* THETA (I). 1*1, NAC 
2.4802E«00 


-1 


NO. OF ACTIVF CONSTRAINTS * 1 OF MHICH 0 ARE VIOLATED 

0 ACTIVE OR VIOLATED CONSTRAINTS ARE SIDE CONSTRAINTS 


ACTIVE CONSTRAINT NUMBERS 

NO. GREATER THAN 6 INDICAIES SIDE CONSTRAINT 


CONSTRAINT PARAMETER. BETA * 2.20652E-11 

OBJ * 2.621 324E*00 NO CHANGE ON OBJ 


BEGIN ITERATION NUMBER 7 
CONSTRAINT VALUES 

-2.3423E+0G 6.7277E-03 -1.9835E+00 -2.623 £-01 -6.4122E-01 -1 

PUSH-OFF FACTORS. THETA (I). 1*1. NAC 
7. 1843E+00 


1 .269 1 f* 00 


•2691E*00 


269l£*00 


FIG. 12 - OPTI MIZ A7ION RESULTS FOR B7AJ1PLS 1 - CCNT 
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NO. OF ACTIVE CONSTRAINTS « I OF WHICH 1 ARE VIOLATED 
0 AC.IVE OR VIOLATEO CONSTRAINTS ARE S.COE CONSTRAINTS 

ACTIVE constraint numbers 

NO, GREATER THAN 6 INDICATES SIDE CONSTRAINT 


CONSTRAINT PARAMETER. BETA » 7.7lS6lE*0i 

S-VECTOR 

2.2189E.01 7.9705E+00 

SLOPE * 7.06315E+01 ALPHA * 2.49714E-04 

OBJ * 2.638962E*0G 

X-VFCTOR 

7.8936E-01 A.0633E-01 


BEGIN ITERATION NUMBER 6 
CONSTRAINT values 

-2. 3333E*00 -4.2505E-11 -1.9775E*00 -2.6666E-0 1 -6.4A19E-01 


PUSH-OFF FACTORS. THETA (I), J«1,NAC 
l.OOOOE.OO 


NO. OF ACTIVE CONSTRAINTS « 1 OF WHICH 0 ARE VIOLATED 

0 ACTIVE OR VIOLATED CONSTRAINTS ARE SIDE CONSTRAINTS 


ACTIVE CONSTRAINT NUMBERS 

NO. GREATER THAN 6 INDICATES SIDE CONSTRAINT 
2 

CONSTRAINT PARAMETER. BETA ■ 3.49062E-IQ 

OBJ « 2.638962£*0O NO CHANGE ON OBJ 


-1.2669E*00 


FIG 


12 - OPTIMISATION RESULTS FOR EXAMPLE 3 - CONT 




final optimization information 

OBJ * 2.638962E*00 

X-VECTOR 

7.0936E-O1 4.0633E-01 

CONSTRAINT VALUES 

-2.3333E*00 -4.2505E-1I -i.9775E*00 -2.6686E-01 -6.AM9E-01 

ACTIVE CONSTRAINT NUMBERS 

NO. GREATER THAN 8 INDICATES SIDE CONSTRAINT 

TERMINATION CRITERION 

ABS(CT) IS LESS THAN CTMIN 


NUMBER OF ITERATIONS * 

8 



OBJECTIVE FUNCTION MAS EVALUATED 

26 

times 

CONSTRAINT FUNCTIONS MERE 

EVALUATED 

26 

times 

GRADIENT OF OBJECTIVE MAS 

CALCULATED 

6 

times 


I .266S>E*00 


FIG. 12 - OPTT^TZATICff RESULTS FOR EXAMPLE 3 - CCNT 
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APPENDIX A 

SDK MARY OF PARAMETERS USED BY C3NMIN 


COMMON BLOCKS: 


/CMhN1/ IPRINT,NDV,ITMAX,NC0N,N3IDE,ICNDI'’ N3CAL NKj’ 

zuAUrUN r LINOBJ r ITFM # ITER r NEAL(4) 

COMMON /CNMN2/ X (N 1 ) , DF { !1 1 ) , G (N2| , ISC ( N 8) , I C { \ 4 ) ’ n-4 \ n 

COMMON /CNMN3/ S (N3) , G 1 <N7) . G2 (N2 [, C ( 55 ). MS 1 <N6 ' 1 n' 

COMMON /CNMN4/ VLB ( N 1 ) , V UB ( N 1 ) , SC AL (N 5) * ' 

CALL STATEMENTS: 


CALL CONflIN (SUB1,CBJ) 
CALL {SUB1 (INFO, OBJ) 


PARAMETERS DEFINED IN THE MAIN PRG3RAM. 

PABAH. DEFAULT DEFINITION 


IPRINT 

HDV 

ITHAX 10 

NCON 

NSIDE 


ICNDIR NDV+1 


NSC AL 


Print control. 

Number of decision variables, x (I) . 

Maximum nuaber of iterations in the 
©ini miration process • 

Nuaber of constraint functions, G(j). 

Side constraint parameter. NSlDL.or.C 
indicates that lower and upper bounds are 
imposed on the decision variables. 
Conjugate direction restart paraaeter. 
Restart with steepest descent move every 
ICNDIR iterations. 

Scaling control parameter. NSCAL.LI. D , 
user supplies scaling vector. NSCAL.Eg.O, 
no scaling. NSCAL.3T.0, automatic linear 
scaling every nVSCAL iterations. 
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PABAM. 

NPDG 

fd:h 

FDCHN 

CT 

CTHIN 

CTL 

CTLHIN 

THETA 

PHI 

MACMX 1 
DEL PUN 


DEFAULT DEFINITION 

Gradient calculation control parameter. 
C.01 ftplati ve change ir. decision variable. 

X(I), in calculating finite difference 
gradients. 

C.01 Kiniaua absolute step in finite difference 

gradient calculations. 

^ Constraint thickness paranetar. 


C. 004 


- 0.01 


0.C01 


1.0 


5.0 


0.001 


Mini.ua absolute value of cr considered m 
optimization process. 

Constraint thickness parameter for linear 
and side constraints. 

Minimum absolute ’alue of CTL considered 
in optimization process. 

Mean value of push-off factor in aethod of 
feasible direct ioas * 

Participation coefficient, used if one or 
more constraints are violated. 

1 pins user's best estimate of the aaxiaum 
number of constraints (including side 
constraints) mhich uill ba active or 
violated at any time in the ninimizatior 
process. 

Minimum relativa change in objective 
function, OBJ, to indicate convergence. 
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PARAH. DEPAULT 


definition 


DA3PUN 0.C01* 
Init. OBJ 

LINOBJ 


ITBH 3 


X(II1) 
VLB (N 1) 

VUB(NI) 

SCAL (MS) 
ISC (N 8) 


nmiauiB absolute change in objective 
function, OBJ, to indicate convergence. 
Linear objective function identifier. 
L 1N0BJ* 1 if OBJ is specifically known to 
be linear in X{I). LIN03J=3 if OBJ is 
non-linear. 

Nunber of consecutive iterations to 
indicate convergence by relative or 
absolute changes, DELFUN or DABFun. 

Vector of decision variables. 

Vector of lower bounds or, decision 
variables. 

Vector of upp er bounds on decision 
variables . 

Vector of scaling parameters. 

Linear constraint identification vector. 
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PAFARETERS 

PARAMETER 

OBJ 

G(N2) 

DP(N1) 

NAC 

A (N4, N3) 
IC(N4) 


PARAMETERS 

PABARETER 

ITER 
NCAL (4) 

S(»3) 

G1CR7) 

G (M2) 

B (M4, N4) 

CCB 9) 


DEFINED IN EXTERNAL ROUTINE SU61 

DEFINITION 

Value of objective function. 

Vector of constraint values. 

Analytic gradient of objective function. 

Number of active and violated constraints 

(G (J) .GE.CT) . 

Rov I contains analytic gradient of the Ith active 
or violated constraint. 

Identifies vhich constraints are active or 
violated. 


DEFINED IN CONMIN AND ASSOCIATED ROUTINES 

DEFINITION 
Iteration number. 

Bookkeeping information. NCAL (I) gives number of 

times INFO=I during optimization process. 
Direction vector. 

Temporary storage of vectors G and x. 

Temporary storage of vector G. 

Used in finding usable-feasible direction. 

Used in finding usable-feasible direction and for 
teaporary storaqe of vector X. 


APPENDIX B 

CONMIN SUBROUTINE DESCRIPTION'S 


P i q « 5 h 


Following is a list of 


the subroutines 


associated with 


CONMIN. if t;»e array dimensions 
currently used, the common blocks in 
changed accordingly. 


are changed from those 


each routine must 


be 


CONMIN - Main optimization routine. 

CNMII01 - Routine to calculate gradient information by finite 
difference. 


CSMN02 - calculate direction of steepest descent, or conjugate 
direction in unconstrained function minimization. 

CNHH03 - Solve one- dimensional search in unconstrained function 
minimization. 


CNMN04 - Find minimum of one-dimensional function by 
interpolation. 


polynocial 


or 


CMNM05 - Determine usable-feasible, 

usable-feasible, direction in constrained 
minimization. 


icodif i od 
function 


CBM NO 6 - solve one-dimensional search for constrained function 

minimization. 
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CNHNC7 - Find zero of one-dimensional function by poi/r.ouiai 
interpolation. 

CNMNC3 - solve special linear programming problem i 
determination of usable-feasible, or -ouified 

usable-feasible direction in constrained function 
uiniaization. 


CNHN09 - Onscale and rescale decision variables before and 
after function evaluation. 
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